Entanglement entropy for odd spheres 



J.S.Dowker 1 

Theory Group, 
School of Physics and Astronomy, 
The University of Manchester, 
Manchester, England 



It is shown, non-rigorously, that the effective action on Z — factored 
odd spheres (lunes) has a vanishing derivative at q = 1. This leaves 
the effective action on the ordinary sphere as (minus) the value of the 
entanglement entropy. Some numbers are given. 



dowker@man.ac.uk 



1. Introduction. 

Finding the entanglement entropy associated with a spatial surface presents 
an interesting computational challenge, e.g. [1-3], that has attracted some recent 
interest, e.g. [4], not least because of the introduction of holographic notions, [5-7]. 
Most attention is given to even space-time dimensions, d, but interest attaches also 
to odd d, [7,8,4], and this is the case I wish to address in this paper for the very 
special situation of spherical surfaces and conformally invariant scalar fields. 

In a previous work, [9], I have analysed the universal logarithmic term that oc- 
curs in the expansion of the entanglement entropy associated with even spheres in 
the special case when the conically deformed Euclidean space-time is the cyclically 
factored sphere, S d /Z q . Then the spatial surface is a (d — 2)-sphere (with vanish- 
ing extrinsic curvature). Conformal transformations can be used to move between 
different manifolds, in particular between spheres, planes and cylinders. 

A universal term is one that is regularisation independent. In even dimensions, 
for a cut-off regularisation, this takes the form of a logarithm, while, for dimensional 
and ^-function regularisation, it shows up as a pole. The 'problem' is that for odd 
dimensions there is nothing analogous (for conformal fields), indeed, in the last 
two regularisations there are no divergences at all and, connectedly, the integrated 
conformal anomaly vanishes. 

In [5] (see equn.(7.11)), it is suggested that, for odd dimensions, the constant 
term in the entropy, i.e. the one independent of any introduced cutoff, is such a 
universal term. Myers and Sinha, [8] , have taken up this suggestion in their attempt 
to produce a c-theorem in higher dimensions and have commented on its validity. If 
there is no divergence, then the entire expression should, presumably, be considered 
universal. 

If boundaries or singularities exist then the relevant C d / 2 coefficient (determin- 
ing the conformal anomaly) is non-zero and a 'proper' universal term might occur. 
This would not be so for the (periodic) factored sphere, S d /Z g (or 'periodic lune'), 
which was the case discussed exclusively in [9] and will occupy me mostly here too. 

In the present paper I wish to look at the entire expression for the entropy. I 
use (^-function regularisation and there is, in general, a pole term, related to the 
integrated conformal anomaly, C d / 2 , and a finite remainder given by a functional 
determinant. I anticipate that the results have interest beyond their entropic bear- 
ing. 

As before, I will obtain the periodic results by combining those for DN-lunes, 
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which are lunes with Dirichlet (D) or Neumann (N) conditions on the boundary. 
Adding these gives the periodic lune, [10,11]. Geometrically, the DN-lune is ob- 
tained by factoring the sphere by the dihedral group, [q], of order 2q. 2 When 
q = 1, the resulting fundamental domain is the DN-hemisphere. 

In [9], I showed that the conformal anomaly on the even dimensional periodic 
lune, S d /Z ? , i.e. the conformal heat-kernel coefficient, C d / 2 (q), had an extremum 
at the ordinary sphere, q = 1. I will show that the same applies to the effective 
action (or functional determinant) for odd dimensions. 

2. The lune 

The (i-lune can be defined inductively by giving its metric in the nested form 
ds 2 d -iune =d6 2 d + sin 2 6 d d8 2 (d _ 1} _ lune , < d i < it , (1) 

which is iterated down to the 1-lune of metric dQ\ with < Q x < (3. The angle Q x is 
referred to as the polar angle and conventionally written <fi. The angle of the lune 
is fi. 

The boundary of the lune comprises two pieces corresponding to (p = and 
cp = (3. The metric (1) shows immediately that these are unit (d — 1) -hemispheres 
because their polar angle, 9 2 , runs only from to 7r. Conditions, typically Dirichlet 
and Neumann, can be applied at the boundary. The boundary parts intersect, with 
a constant dihedral angle of (3, in a (d — 2) -sphere, of unit radius, which constitutes 
a set of points fixed under 0(2) rotations parametrised by (p. 

It can be seen that the 2-lune submanifold, with coordinates 9 1 and 0, has a 
wedge singularity at its north and south poles of. These poles are at 6> 2 = and 
6> 2 = iv and are the O-hemispheres of a 0-sphere. In the <i-lune, the submanifolds, 
6> 2 = and 6> 2 = 7r, are the (d — 2) -hemispheres of the (d — 2)-sphere of fixed points 
just mentioned. 

All this is for arbitrary angle, (3. If /3 = ir/q, q e Z, the lune can be identi- 
fied with the fundamental domain of the dihedral group action on the sphere. As 
mentioned previously, restricting to just the rotational subgroup doubles the fun- 
damental domain to a periodic lune of angle 2ir/q, the edges now being identified. 
This is equivalent to adding the sets of and D modes, [11] and produces a conical 
singularity of deficit 2tt — 2/3. It corresponds to the introduction of the uniformising 

2 This is generated by reflections in two hyperplanes in the embedding space. 
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angle, 4>' = 2n(j)/2/3, which runs from to 2ir as 4> runs from to 2(3. The defor- 
mation of the (i-sphere, (1), can be traced to that of the 1-lune, or ^/-circle, which 
has unit radius but circumference, 2/3. 

It is sometimes convenient, for visualisation purposes, to think of the (i-sphere 
as embedded in M. d+1 with the lune pictured as the curved, outer surface of a 
'hyper-wedge'. The edges of the lune are great unit (d — l)-hemispheres, half the 
intersections of the sphere with two <i-flats (hyperplanes) mutually intersecting in 
a (d — l)-flat, the 'axis' of the 'rotation' that takes one <i-flat into the other with 
rotation angle, (3. This axis intersects the sphere in two 'poles', isometric in total 
to a (d — 2)-sphere and a singular region. It is a fixed point set for those elements 
of 0(d + 1) that maintain the axis of rotation, i.e. an 0(2) subgroup. 

If j3 = ir/q, q G Z, the lune, and its q reflections in the hyperplanes, tes- 
sellate the (i-sphere. Combining a lune with its neighbouring reflection gives the 
fundamental domain for the rotational (cyclic) subgroup, Z , of the dihedral group. 

To make the geometry a little clearer, I look at the d = 3 case and write out 
the Cartesian embedding of the unit 3-lune, in my notation, 

x 1 = cos 6> 3 

x 2 = sin 6> 3 cos 6> 2 

x 3 = sin 6> 3 sin 6> 2 cos <j> 

x 4 = sin 6*3 sin 6> 2 sin 

with < <p < (3. Translations in <f> correspond to rotations in the x 3 -x A plane with 
the fixed point x 3 = x 4 = giving the set of fixed points on £| as the unit circle, 
x i + x 2 = 1 f° r au openings, (3. 

The example of (3 = n is easiest to picture. The boundary pieces of at 
<p = and = 7r, form, respectively, the x 3 > and x 3 < hemispheres of the 
2-sphere, x 4 = 0. The boundaries of these pieces coincide with the fixed point unit 
circle x\ + x\ = \ and join, with a dihedral angle of n, to form an equator of the 
2-sphere with poles, x 3 = ±1. 

3. Conformal Transformations 

The conformal transformations relevant here are those considered in [12], i.e. 
those between the original sphere, S d , the Euclidean plane, M. d , and the cylinder, Ix 
S d_1 . The latter would give the entropy on a (d — l)-sphere where the separating 
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surface is the equatorial (d — 2)-sphere. This has vanishing extrinsic curvature, 
whereas the second mapping yields the entropy on a flat (d— 1) space with surface 
an ordinary (d — 2)-sphere which does have extrinsic curvature 'generated' by the 
conformal transformation. Actually, this geometrical circumstance was behind the 
calculations of [13,14] and [6] and I am, in reality, going over part of this ground 
again. 

In the particular example of the 3-lune, consider the section i 4 = to be a 
space section. Project it from the S 3 pole x 3 = 1 onto the equatorial plane, x 3 = 0. 
The fixed point circle projects to a unit circle, centre the origin, the inside of which 
is the projection of the hemisphere with negative x 3 and the outside is the projection 
of that with positive x 3 . 

The log coefficient is conformally invariant, and so it is immaterial where it 
is evaluated. However the entropy, being determined by the effective action, will 
generally change (in a well defined fashion). 

The standard anomaly equation for the conformal variation of the effective 
action, W, is, in <i-dimensions, 



Equation (2) can be integrated if the local coefficient is known. For arbitrary 
dimension we do not have this luxury and, in any case, the evaluation is complicated. 
It has been done only for d = 2,3,4 and 6. However, life is made easy in odd 
dimensions, if no boundary is present (as is the case for the periodic lune) for then 
the right-hand side of (2) is sero. Therefore the effective actions on conformally 
related spaces are the same, and likewise for the entropies according to the standard 
recipe given in the next section. 

4. Entropy 

According to the general prescription of Callan and Wilczek, [15], the entropy 
is given in terms of the effective action, W^ B \ on the Euclidean manifold deformed 
by a conical singularity of angle 2n / B as, 




/ 




S = -(Bd B + l)wW 



(3) 



B=l 
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In the present case, the deformed manifold, M. q , is the periodic lune, S d /Z g , 
so that B = q. The field theory under consideration is the scalar field conformally 
coupled in d dimensions and the effective action for this has been calculated in 
[16,17]. The evaluation of the derivative in (3) involves a continuation in q from the 
integers 3 . In our previous work this was straightforward as the relevant quantities, 
the heat-kernel coefficients, were polynomial in q. For the effective action, the 
matter is a little less obvious. 

In terms of the (^-function, £(s 7 g) ? on JA , the unrenormalised effective action 
is, generally, 

W M = I lim - V(0, q) + C(0, q) log L + X , (4) 

Z s— >-0 S z 

where X is a possible finite correction and £(0, q), = C d / 2 (q), the conformal anomaly. 
The third term, where L is a scaling length, can be considered a concomitant of 
the ultraviolet pole divergence, (in this version of (function regularisation) . Taken 
together they correspond to the log term that arises in the cut-off method. Conven- 
tionally, the pole divergence would be removed by renormalisation. The 'area law' 
for the entanglement entropy is not recovered nor the other divergences encountered 
in a cut-off or a lattice approach. This is not a worry as these are non-universal 
terms anyway. 

Generally, when ((0) vanishes, its role is, in some ways, taken over by the 
derivative, C'(0), this then also being conformally invariant. Hence, as a working 
hypothesis, when treating the periodic lune, I set X to zero (it is a conformal 
invariant) and take the entropy, 



(5) 

q = l 



to be universal, in accordance with [5] and [8] and now have to address the calcu- 
lation of the derivative. 

5. The log/divergence term 

As mentioned, in the present instance where d is odd, the conformal anomaly 
on the periodic lune is zero, 

C d/2 (q) = C? /2 (q) + C? /2 (q) = , 



3 The theory for any real q is easily worked out using eigenfunctions. 
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and there is no pole nor log term. Nevertheless, for completeness, I examine the 
individual coefficients as they would be required for analysis on a DN-lune singly. 

The required C d / 2 (q) coefficients are given by the expressions in [9], where the 
original references can be found. One has, 

C?/M = -C° /2 (q) = ^ (Bf\d/2 - 1| q, 1) + Bf\d/2\ q, 1)) , (6) 

The symmetry of the generalised Bernoulli polynomials, 

S J0(d _ 1 + q - x \ q , 1) = (_!)- Bi d \x\ q, 1) , 

should be noted. 

For odd d the right-hand side of (6) simplifies to a constant, independent of q 
and, therefore, 

= -cS,M = s J°( d / 2 - 1) , (7) 

which is proved in the Appendix. 

Thus again, as in [9] , the g-derivative of the conformal anomaly is zero (trivially 
in this case) at q = 1. Numerical values are 



= ± ld^ B i B d-i\d/2-l) (8) 

1_ _L9_ 17 271 _ 367 3233 
^~ 48' ~ 1440' 11520' 120960' ~ 1935360' ~ 7257600' " ' 

for d = 3, 5, . . .. 

This (/-independence can be seen in a different way. For a <i-dimensional man- 
ifold, M., the heat-kernel coefficients, C k , in general take contributions from all 
submanifolds, of codimension down to d. It is easily checked that the integrand 
of the half-integral coefficient, C fc / 2 ? f° r even codimension, has dimensions of an 
inverse odd power of length. Its construction must then necessarily involve powers 
of the extrinsic curvatures but these all vanish for the fundamental domains of the 
group action on the sphere, which applies here. In particular, the codimension two 
contribution to C d / 2 is zero and this is the part that involves the 'dihedral' angle, 
7r/q between the boundary parts of the manifold, M.. 

For example, in three dimensions, C 3 / 2 then has just an area contribution 
and a trihedral corner, or vertex, contribution, [12], but, in the case of a dihedral 
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action, the vertex degenerates into a dihedral one, and its contribution vanishes. 
The explicit local form of the area part remaining is, 



Coin — ± 



3/2 384tt 

where R = 2 is the intrinsic curvature of the boundary part dM i: a 2-hemisphere. 
These parts add to a full S 2, s worth of boundary (for all q) and (9) evaluates to 
±1/48, in agreement with (8). 

6. The derivative term 

I now consider the universal term in the entropy (5) (actually the whole en- 
tropy). For the periodic lune, the (-function is the sum of the N and D (-functions, 

C(s»?) = (n( s ,q) + (d( s ,q) ■ 

The dihedral (-functions have been calculated in [16], and in [17], in terms of 
the Barnes (-function, ( d (s, a\ d). In particular the derivatives at zero are, 4 



&(0, q) = Cd(0> d / 2 \ V, 1) + Cd(0> d/2-l\ q, 1) + M(a N , q) 

Cd(0, q) = Cd(0, d/2 + q\ q, 1) + ( d (0, d/2 - 1 + q\ q, 1) + M(a D , q) , 

where M(a, q) is given by, 



(10) 



(d-l)/2 H Q 

M(a,q) = - £ y^r N 2r(d,q), (11) 



with the odd harmonic number, 

1 



The N L are the residues of the Barnes function, 



N, 

Q(s + l,a\d)-> — + R t , ass^O. (12) 
s 



4 The first two terms on the right-hand side are suggested by the factorisation of the eigenvalues 
with the final term indicating a multiplicative anomaly. 
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They are given by generalised Bernoulli polynomials (see later) and depend on d, 
the argument, a, and the parameters, d, which here reduce, in effect, to just the 
number, q. The dihedral Neumann and Dirichlet arguments, a N and a D , are given 
by (d — l)/2 and (d — l)/2 + q, respectively. 

The two lines in (10) give the derivatives of the ^-functions for Neumann and 
Dirichlet boundary conditions on the edges of a lune of angle n/q. As described 
earlier, adding them gives the ^-function for a (periodic) lune, or cone, of angle 



2n/q that is, for S d /Z g which is the object of most interest. 

According to (3), I need the derivative of (10) with respect to q. This is non- 
controversial for the M terms as they are rational in q. For the Barnes function, I 
will assume that the derivatives with respect to s and q commute and that q can be 
continued to 1. This can be justified from the contour expressions for the Barnes 
multiple functions (see the Appendix). 

There now follows a technical discussion of the derivatives of the Barnes func- 
tion needed in (10). 

7. Derivatives of the Barnes ^—function 

Starting from the formal definition 5 , 

oo ^ 

Q(s, a\ q, 1) = V — — (13) 

^ {a + qn 1 + ...n d ) s 

simple manipulation leads to the derivative at q = 1, 
d q Q(s,a\q,l 



^ = -- C d (s, a\ 1) + s- Us + 1, a\ 1) , (14) 

assuming that a is independent of q. Relevant for this is the important recursion, 

( d (s,a + q\q,l) = ( d (s, a\ q, 1) + Cd-i(s, a\ 1) , (15) 

which, looking at (10), relates the N and D expressions. In particular it means that 
the Dirichlet contribution to the g-derivative equals the Neumann one. 

The derivative with respect to s and the limit s — > can now be taken in (14) 
with a view to substitution into (10). Thus, interchanging the implied limits, 

q =i a a (16) 



9„ft(0,a| 9 ,l) 



5 The continuation in q is visually obvious. 
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R 1 is the remainder at the s = 1 pole, (12), when all parameters equal 1. Henceforth, 
I do not display dependence on these parameters unless necessary. 

Before putting things together, I give some details on the residue and remain- 
der. The latter is given by Barnes as 

R r (d, a) = (-iy ((73T)T ^ (a) - N rV> q) H *-i) ( 17 ) 
where H r is the usual harmonic number, 

r 1 

H r = T » #o = . 



k=i 



and the ^-functions are defined in terms of the multiple T-function, 

#(«) = ^logr d (a) 



(18) 



Hence the residue R^d, a) = —ip^\a) = —ip d (a) and the derivative (16) becomes 



dq(' d (0,a\q,l) 



q=l 



dd , Bf ( a)--^M«)- 



(19) 



It is shown in the Appendix that this can be rewritten more compactly as 



d q C d (0,a\q,l) 



= ip d+1 {a + l) . 



5=1 



(20) 



This result will be employed shortly. 

I emphasise that I use Barnes' definitions of the multiple functions. 



8. Derivative of the multiplicative anomaly term 

I now compute the derivative of the multiplicative anomaly, M, terms in (10). 
I do not need to do this for the periodic lune, as these terms cancel on addition, as 
shown next. However I give the results, again for completeness. 

I need the residues, , given in terms of generalised Bernoulli polynomials as 



1 1 n {d) 



(21) 



(l-l)\(d-l)\q 



-B^d-l + q-alq,!)- 
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If I = 2r, as in (11), and d is odd, the minus sign in (21) shows that the M 
terms in (10) are opposite for N and D conditions and would cancel if added, as 
advertised and as noted in [16-18]. 

I show that the (/-derivative of M(a N ,q), (11), vanishes at q = 1 which cir- 
cumstance follows directly from the general result, [19], 

Qd q ~ q BflM q, 1) = -1 B%?\a + q\ q, q, 1) . (22) 

For, setting a = a N = (d — l)/2 and q= 1, the right-hand side equals 

-B<££\{d +l)/2) 

which vanishes from a well known property of the Bernoulli polynomials since d — 2r 
is odd. 

It is interesting, and potentially more significant, to note that the multiplicative 
anomaly, M(a N , q), is, in fact, independent of q, if d is odd, just like the conformal 
anomaly. For d = 3, 5, . . ., one has, 

1 5 2303 142601 

{a N ,q)---, — , - 691200 ' 203212800''"' ^ ' 



9. The q— derivative on shell 

Formally differentiating the full ^-functions, (10), with respect to q at 1, using 
(20) and the result of the previous section, I find for the total Neumann ^-function, 
(the Dirichlet value is the same) that, 



= Vd+i(d/2) + V*fi(d/2 + l). (24) 

q = l 



d q C' N (0,q) 

The expression, 

M°) = ^9^-i W - (^yy E ^ S^iCd - a) sfV) , (25) 

for the multiple ^function in terms of the standard digamma function, ifj(a), is 
derived in [20], where other references can be found. 

Then, taking into account the antisymmetry for odd d, 

B { f +1) (d/2) + B ( f +1) (d/2 + 1) = , (26) 

10 



and the standard formula for the digamma function, 



V>(l/2-n) = V(l/2 + n) = -7-21og2 + 2#°_ 1 , tiGZ, 

the transcendentals 7 and log 2 cancel leaving the purely algebraic expression, 

= -^f +1 \d/2 + l)- 

, u d^ ( ^ A k Md/2)B^- k) (d/2 + l) 

f^y i} dTT^i (27) 

k=l 

1 J J d+l-Jfe 

fc=l 

which, unsurprisingly perhaps, evaluates to zero, by machine, dimension by dimen- 
sion. 

An expression for the second (/-derivative of the effective action is derived in 
the Appendix. Its numerical evaluation shows an alternating sign, being positive 
for d = 3 (corresponding to a minimum for £'(0, q) at q = 1). 

10. The entropy 

According to this result and the definition, (5), the entanglement entropy re- 
duces to just half the quantity, C'(0, (?) | x , obtained by adding the ordinary ND- 

hemisphere effective actions. These have been given above, (10), and can be formally 
expressed in terms of multiple T-functions, 

CD{0 , 1) = loe EMl±Wl. M(aN , lh (28) 

Pd 

where p d = r d+1 (l) is the multiple modular form, [19], and the M term is given by 

(23). 

The full sphere effective action has a longish history, some of which is detailed 
in [18]. From the sum of the N and D expressions, (28), (see [17] equns.(22),(23)), 
I find, 

= r d (d/2-i) r<l (d/2 + i)rg(d/2) 

Pd 
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Explicit formulae can be obtained by expanding the Barnes (^-function in terms 
of the Hurwitz ^-function with coefficients related to Stirling numbers, [21]. I just 
give some samples on the full sphere, 



C'(0, 1) = -| Cfl(-2) - \ log 2 » -0.127614 , d = 3 

= ~ Cfl(-4) - ^ Cfl(-2) + ^ log 2 « 0.011486 , d = 5 . 
Another form follows by rewriting (29 as 

r d+1 (d/2-i)r d+1 (d/2) 



(30) 



C'(o, 1) = log 



1 d+1 

(d/2 + 2)r d+1 (d/2 + l) 

d+1 rl {d-l)/2 (31) 



2(— l) d+i f 

= -i— ^ J dz Tvz tanyrz JJ (^ 2 - (j - 1/2) 2 ) 

which is derived in [20]. It corresponds to the k = 1 case of the GJMS operator, 
P 2 fc? which is just the usual conformal Laplacian. The integrand can be expanded 
and contact made with (30) but for numerical purposes (31) is adequate and I find 
for C'(0,1), 

-0.001595, d = 7 

0.000262 , d = 9 (32) 

-0.000047, d=ll. 
The values alternate in sign, like the conformal anomaly in even dimensions. 



11. Conclusion 

The main technical result in the preceding is the vanishing of the g-derivative 
of the effective action on the odd dimensional orbifold S d /Z q at q = 1, showing that 
the entanglement entropy associated with a hyperspherical (d — 2)-submanifold is 
essentially just the effective action on the ordinary d-sphere, on certain assumptions 
regarding universality. 

I have not yet been able to prove the responsible Barnes-Bernoulli identity, 
which, in terms of multiple functions reads, 

V> d+1 (d/2) + V> d+1 (<V2 + l) = 

for all odd d. 

The entropic significance of the numerical values on shell, (30), (32), is unclear 
to me. Perhaps more interesting would be results off criticality obtained, say, by 
adding a mass term. Extension to other fields is also indicated and also to other 
symmetric spaces. 
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Appendix 

In this appendix I give some relevant material concerning the Barnes functions 
and start with the proof of the (/-independence of the N and D C d / 2 coefficients, (7), 
which I will do algebraically. I seek to show that the bracket in (6) is linear in q 
and write down the standard polynomial, 

Bl n) (*|«, 1) = Q B s BtS 1] (x) , (33) 

s=0 

noting that the only term odd in q comes from the s = 1 term in the sum. Next I 
apply the formula ([22] p. 167), 

Bi n \x + q\q,l) = Bi n \x\ -q,l) (34) 

to get 

B f\d/2 — 1] — q,l) = B [ d d \d/2 - 1 + q\q, 1) 

= -Bf (d/2|g,l) 

after using symmetry for odd d. Hence the bracket in (6) is twice the odd part of 
B d d \d/2\ q, 1) obtained from (33), thus proving (7). 

Contour integrals provide an alternative starting point. Barnes [19] p. 406 gives 
an expression for the derivative at zero, 6 

e" az (log(-z) + ^)dz 



CM a\q,q,l) = ^j c z(1 _ e -,, )(1 _ e - q *, )(1 _ e -,)d-2 ^ 



and so 



n ,, i . if e-^ a+ ^ z (\og(-z)+^)dz , 
d Cd 0, a\q, q', 1 = — / - V , f- ; . 36 

To compare, Barnes also gives the integral for the -i/^-function, obtained simply 
by differentiating (35) with respect to a, 

/ / I if e~ az (log(- z) +-f)dz 



6 „/ 



q is introduced for notational flexibility. 
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so we obtain very simply (set q' = 1 in (36) and q' = q in (37)), 

dqC'd(0, a\q, 1) = V>d+i( a + q\ q, q, 1) 



(38) 



Evaluating at q = 1 yields, 



d q ('(0,a\q,l) 



q = l 



\d-i 



a . , . (-I)" - D (d) 



(39) 



after using Barnes' recursion, 



V'd+iCa + 1) = - ^ C d (0, a) , 



(40) 



and so I have regained the expression (19), derived earlier by a different method. 7 
Higher derivatives can be deduced. For example 



d 2 q (' d (0, a\q, 1) = 2^ d %(a + q\ q, q, q, 1) + ^(a + 2q\ q, q, 1) 
whence, setting q = 1, 



(41) 



^(0,^,1) 



_ o„/,(2) 



2V>£ 2 (a+l) + v2!i(<» + 2) 



5=1 



da 2 



(42) 



(2^ +2 (a+l) + V d+ i(a + 2)) 



expressed via the trigamma function. This process can be continued easily. 

Algebraically, the right hand side of (24) vanishes because of the factorisation 



V> d+ i(a) + + 1) = (a - d/2) (p d {a) ^(a - d) + 



Qd(a) 
a — d 



for odd d where P d (a) and Qd( a ) are polynomials such that there is no pole at 
a = d. 



7 One could turn this around and use this development to prove the recursion, (40). 
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